Vortex Contribution to Specific Heat of Dirty d-Wave Superconductors: Breakdown of 

Scaling 
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We consider the problem of the vortex contribution to thermal properties of dirty d-wave supercon- 
ductors. In the clean limit, the main contribution to the density of states in a d-wave superconductor 
arises from extended quasiparticle states which may be treated semiclassically, giving rise to a spe- 
cific heat contribution SC(H) ~ H 1 ^ 2 . We show that the extended states continue to dominate 

the dirty limit, but lead to a HlogH behavior at the lowest fields, H c i^H <C H c 2- This crossover 
may explain recent discrepancies in specific heat measurements at low temperatures and fields in 
the cuprate superconductors. We discuss the range of validity of recent predictions of scaling with 
H 1 ' 2 /T in real samples. 



PACS Numbers: 74.25.Fy, 74.72.-h,74.25.Jb 

Introduction. With the growing consensus that the 
symmetry of the cuprate superconductors is d-wave 
has come a renewed interest in the properties of the vor- 
tex state in "unconventional" superconductors with or- 
der parameter nodes. Many of the basic ideas about how 
this state differs from the conventional Abrikosov state in 
classic superconductors were worked out already in the 
context of rotating 3 He and heavy fermion superconduc- 
tors. Recently, however, a number of novel features of the 
problem peculiar to those systems with Dirac spectrum 
(line nodes in 3D or point nodes in 2D with order param- 
eter vanishing linearly with angle on the Fermi surface) 
have been pointed out. Volovik showed that, in con- 
trast to conventional superconductors, extended quasi- 
particle states with momentum k near order parameter 
nodal directions k„ dominate the density of states at 
zero energy. This leads to a specific heat which varies 
as SC(H) ~ H 1 / 2 , in contrast to classic superconduc- 
tors, where localized quasiparticle states in vortex cores 
lead to a scaling of 8C(H) ~ H since the number of vor- 
tices scales proportionally to the field. Simon and Lee 
H then showed that thermal and transport properties 
exhibit a scaling with H 1 ^ 2 /T, again arising simply from 
the low-energy Dirac form of the electronic spectrum. 

The predicted proportionality of the electronic specific 
heat to \[H was in fact identified in measurements on 
high quality single crystals by Moler et al. , Q one of the 
crucial early experiments lending credence to the d-wave 
hypothesis. However, the interpretation of the observed 
s[H dependence has been questioned by Ramirez [f| who 
points out that there are well-known cases where classic 
superconductors show a similar "nonanalytic" behavior 
sufficiently close to the lower critical field H c \. Further- 
more, experimental results of Fisher et al |(| and Revaz 
et al. cannot be well fit by a VH form. 

The above scaling predictions hold, strictly speak- 
ing, for clean d-wave superconductors and for energy 
scales small compared to the maximum gap scale Ao. 
To make realistic predictions for experiments, deviations 
from scaling due to disorder and other real-materials ef- 
fects must be accounted for. In this work, we study the 



effects of disorder, primarily in the unitarity scattering 
limit thought to be relevant to the cuprates, |^|| and ask 
how the scaling predictions of refs. ||,|| break down. The 
treatment is similar in spirit to the crossover phenomena 
studied in the context of the nonlinear Meissner effect 
in d-wave superconductors by Yip and Sauls fl(i[| In the 

present work we are concerned with fields H~H c i , how- 
ever, and discuss, in a crude way, the influence of the 
structure of the vortex state itself on thermodynamic 
bulk measurements. We analyze existing experiments 
and point out how they may be reconciled with the d- 
wave hypothesis and the ideas of Volovik by properly 
accounting for the effects of disorder. 

Treatment of extended quasiparticle states. In an ex- 
ternal field, one should in principle solve the Bogoliubov- 
de Gennes equations or equivalent for the fully self- 
consistent, spatially dependent mean fields and quasi- 
particle amplitudes. Imposing a vortex- type boundary 
condition on the order parameter around an isolated 
singlularity leads to a uniform phase winding with cor- 
responding superfluid velocity v s = (h/2mr)9, with 9 
the azimuthal angle in real space. The order parame- 
ter magnitude Afc(R) is supressed near the vortex core 
over a length scale of order the coherence length £o, an d 
has a fourfold symmetry in real space. [ pd| The quasi- 
particle wavefunctions and quasiparticle density of states 
have also been claimed Q to display fourfold symmetry. 
These excitations are of two types: bound states local- 
ized in the vortex cores, and extended states which evolve 
smoothly into the bulk zero-field quasiparticle states at 
large distances from the vortex. In a classic superconduc- 
tor, the low-temperature entropy is dominated by the 
core bound states, since the extended states are fully 
gapped and therefore essentially depopulated. The core 
levels are furthermore separated by a typical spacing 
A 2 /Ep, where Ao is the gap maximum and Ep is the 
Fermi energy. This can be rather large in short-coherence 
length superconductors like the cuprates, such that only 
one or a few states may actually be bound. Even if 
treated as a quasicontinuum, in a d-wave superconduc- 
tor the bound states may be shown to contribute less to 
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the entropy than the extended states with momenta near 
the bulk gap nodal directions, by an amount which 
diverges logarithmically with the vortex size or intervor- 
tex separation. Since nonmagnetic disorder is pairbreak- 
ing in a d-wave superconductor, the contribution of low- 
energy extended states to the entropy will be still greater 
in the dirty systems we consider. 

For the above reasons, it appears in the d-wave case to 
be a good approximation to ignore the core excitations 
in the calculation of bulk properties. The extended exci- 
tations are treated here by a method originally proposed 
by Maki and Tsuneto jl3| for classic superconductors and 
applied to the d-wave nonlinear Meissner effect by Yip 
and Sauls. jlO| The semiclassical approximation treats 
the quasiparticle states as plane waves of energy Dopplcr 
shifted by u> — > ui — v s ■ k. The single-particle matrix 
Green's function for the pure system is therefore 



ff (0) (k^;v s ) 



(w - v s • k)r + A k T! + £ k r 3 
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where the 7V are Pauli matrices in particle-hole space. In 
(1), £fc is the usual single-particle band measured rela- 
tive to the Fermi level, and A k = Ao cos 2<fi is the bulk 
d x 2_ y 2 order parameter over a cylindrical Fermi surface, 
taken independent of position in real space. Thus, we 
neglect both the order parameter suppression and the 
spatial variation of the quasiparticle amplitudes near the 
vortex core over the coherence length £ = hvp/irA . 
This is justified provided we confine our interest to fields 

H such that H c1 ^Sh <C H c2 : the coherence length £ 
will then be much smaller than the other relevant length 
scales, specifically R, the intervortex distance and A, the 
penetration depth, in a strongly type-II system. 

Disorder is introduced in the self-consistent t-matrix 
approximation |l4|,li1| through the averaged self-energy 
T,q(uj) = TGqTo I Jcr — Go) , where T — rii/nNo is an impu- 
rity scattering rate proportional to the concentration rii 
of point potential scatterers, c = cot 5q is the cotangent of 
the s-wave scattering phase shift Sq , and No is the density 
of states at the Fermi level. We focus here on the strong 
scattering limit So ~ 7r/2, but results are easily obtained 
for general phase shifts. The averaged integrated Green's 
function is Gq(uj ) = (7rA r ) _1 S k ^ Tr r g(k, uj), which for 
a simple d x 2_ y 2 state leads to Go(ui) = — i(2/7r)K(A /(D). 
This form of the self-energy leads via the Dyson equation 
to an averaged propagator g(k, uj) identical in form to (1) 
but with ui replaced by Co — uj — So(w). Note the Green's 
function must be determined self-consistently, and de- 
pends on the single energy variable uj— v s k. It is further- 
more important to observe that the most general form of 
the propagator would include renormalizations of both 
the single-particle energy £ k and the order parameter A k 
as well. Corrections due to the former vanish identically 
for one-particle properties like the density of states if the 
system is particle-hole symmetric. jl6| Corrections due to 
the latter, which vanish for a d x 2_ y 2 order parameter in 
the H = case, are nonzero in general in finite field, due 



to the dependence of the shift v s • k on the angle <f> over 
the Fermi surface. However, for low energies such that 
only quasiparticles in the neighborhood of the nodes are 
relevant, this renormalization may be shown to be small, 
and we have neglected it here. 

Density of states at zero energy in nonzero field. We 
first present results for the magnetic field dependence of 
the density of states at zero energy, N(0;H). This is 
an experimentally accessible quantity, as it scales with 
the linear-T term in the low-temperature specific heat in 
the superconducting state, 7 e ; = 7r 2 iV(0)/3. The ori- 
gin of this term has been controversial, and may re- 
sult in part from nonelectronic two-level systems away 
from the planes. The zero-field residual density of states 
used in this work arises solely through disorder in the 
electronic system, and must be assumed to represent a 
lower bound to the true residual density of states. To 
find N(0;H), we average the propagator over a vortex 
unit cell, iV(0; H)/Nq = (— ImGo(o;; v s ))h, where for any 
/(v.) we define (/(v s ))„ = A~\H) J ccll d 2 r /(v.). For 
simplicity, we take the cell to be circular, of area A ~ 
ttR 2 . Here R = ^o(7r/2) 1 / 2 a- 1 {H c2 /H) 1 / 2 is the inter- 
vortex spacing, and a is a constant of order unity depen- 
dent only on the vortex lattice geometry. Note that in- 
creasing the magnetic field does not affect v s , but merely 
decreases R. Effects of the actual spatial dependence of 
v s in the lattice phase can be easily incorporated in the 

theory. In the clean limi t and for H c ilSH <C H c2 , we 
obtain N{0;H)/N o ~ y/8/^a(H/H c2 ) 1/2 . This is essen- 
tially the result obtained by Volovik, || who linearized 
the gap and the electronic spectrum around the nodes, 
and may be derived easily by recalling that for the d x 2_ y 2 
state N(lo)/Nq ~ uj/Aq for u> <C Ao, replacing u by 
v s • k, and performing the spatial integral over the cell. 
This approximation must therefore fail at low energies, 
where the density of states in the disordered H = sys- 
tem has the form N(uj) - iV(0) + buj 2 , with AT(0) related 
to the zero-energy quasiparticle scattering rate 270 by 
2V(0) = (27o/7rA )log(4A /7o). In the unitarity limit, 
70 is well approximated by 70 ~ 0.61-v/T Ao for small 
concentrations. 

As in the zero-field case above, the residual density of 
states depends on the quasiparticle lifetime, which be- 
comes, however, a local quantity in the presence of the 
superflow field v s (r). At low energies, substitution of the 
form (I>(r) = ^ + 27(1") into So(^) yields for v s kp, 7 <C Ao 
but arbitrary v s kp/j, 



-^ = -i fin 



4A 



2 A 



+ 



■ tan 



v /7 2 + (v,-k„) 2 ^ 
'v., • k,. 
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where k„ is the nodal direction. In the "dirty limit", 
(H/ H^) 1 ! 2 Ao < 70 < Ao, the spatial integrations in 
(— ImGo)ir can be performed, yielding 
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FIG. 1. Normalized density of states N(0;H)/N for 
T/Ao = 0.1, 0.01, and 0.001 in unitarity limit c = 0. Nodal 
approximation (solid lines), exact result (crosses), clean, low 
energy approximation (dotted line). Data from Fisher et al. 
H (circles); Moler et al. [Q] (untwinned sample, squares), as- 
suming # c2 /a 2 =300T, 7n = 15 mJ-mol-K 2 . 

To illustrate the deviations from the clean limit, low- 
energy result, we plot in Figure 1 a full numerical evalua- 
tion of the density of states N(0;H). For the cleanest case 
plotted, the predicted scaling is attained over a substan- 
tial range. Deviations from the \f~fi behavior occur below 
a lower crossover field scale H* = (j /aA ) 2 H C 2, corre- 
sponding to the smearing of the linear density of states 
due to scattering by impurities. We note that excellent 
approximate results may be obtained, with a substantial 
increase in speed of numerical evaluation, by replacing 
v s ■ k in g(k, uj;~v s ) by its value near the nodes, v s • k„, 
so that G (uj; v s ) ~ J2a=± K ( A o/(^ - «v s • k n )). 

We use this approximate form, whose validity is shown 
in Figure 1, in all further calculations. 

In Figure 1 we have also shown experimental data 
of Moler et al. [Q and Fisher et al. We have plot- 
ted published data for the linear electronic specific heat 
coefficient 7e ; without attempting to subtract residual 
"extrinsic" contributions possibly due to 2-level systems, 
etc. We observe that the data for the untwinned YBCO 
crystal of Moler et al, which are claimed to follow the 
clean-limit prediction of Volovik, Q are actually consis- 
tent with a slightly dirty d-wave superconductor as well if 
the entire residual contribution to iV(0; H) is attributed 
to impurities (analysis of NMR and penetration depth 
data on similar samples makes this interpretation un- 
likely, however. Q). No good fit to the clean limit pre- 
diction was obtained for the sample studied by Fisher 
et al. From the Figure it is clear that this discrepancy 
might be explained by simply assuming the Fisher et al. 
sample contains roughly an order of magnitude more de- 
fects than the Moler et al. sample. Clearly systematic 



disorder studies would be very useful in resolving these 
questions. 

Density of states at finite frequency. Kopnin and 
Volovik JTt| have observed that the field-dependent part 
of the density of states, 8N{uj\ H) will diverge at low fre- 
quencies as 1 /oj down to a lower cutoff of order the av- 
erage quasiparticle energy shift Eh = a(H / H^) 1 / 2 Aq. 
In the current formalism, this result is recovered by not- 
ing that the field-dependent part of the local density of 
states is given by (— Im[Go(w; v s ) — Go(co)])h, which for 
7o < E Hl lo <C A yields 




o — 
q I 1 1 1 

°0.00 0.05 0.10 0.15 



w/Ao 

FIG. 2. Density of states at a(H/H c2 ) 1/2 = 0.025 vs. freq. 
oj/A for r/T c =0. ,0.0001,0.001, 0.01 (solid lines). Interme- 
diate-frequency asymptotic result 8N ~ a 2 TvAoH/(4:UjH C 2) 
(dashed line). 

obtained in the limit of large frequencies, SN(uj; H)/Nq ~ 
a 2 irAoH/(4:LjH C 2) . At energies smaller than the average 
quasiparticle energy shift Eh, this divergence is cut off 
asSN(u;H) ~ N(0; H)(l-irx/i)). Both these limits arc 
clearly visible in the clean case shown in Figure 2. The 
1/uj divergence may also be cut off in a different fash- 
ion if the impurity scale 70 exceeds the magnetic energy, 
Eh < ijJ < 70, as also shown (the maximum in each curve 
corresponds roughly to the scale 70). In dirty samples, 
such that 70 becomes an appreciable fraction of the gap 
scale, the l/to behavior will be unobservable. 

H 1 / 2 /T scaling of specific heat. The specific heat is 
now easy to calculate by differentiating the entropy of a 
free Fermi gas of quasiparticles with disorder- and field- 
averaged density of states N(uj; H); one finds at low tem- 
peratures: 
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c 



duj 



-df 



dui 



N(u;H) 



(5) 



N(0;H)^-T 
T 2 



T < max[j 0> E H ] < A 
7o,£ ff « T « A 



A T 2 term characteristic of the pure d-wave system in 
zero field is present whenever both the impurity and mag- 
netic field scales are smaller than the temperature. If one 
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FIG. 3. Electronic specific heat C e i(H)/(-y n T c ) at 
T = 0.1T C normalized to normal state value, 7„T C vs. square 
root of field, a(H/H c2 ) 1/2 for T/T c = 0.1, 0.01 and 0.001 (solid 
lines). Clean limit (small filled circles). Asymptotic low-T 
clean limit (dashed line). Zero-field linear term N(0; 0)T 
(filled triangles) for values of T shown. 

plots the specific heat at fixed temperature, as in Fig. 3, 
this results at low fields in deviations from y/~H behavior 
due both to impurity induced residual density of states 
and, in cleaner samples, the T 2 term. The latter effect is 
the origin of the saturation of the progressively cleaner 
curves in the figure to a non-Vi? behavior. 

Finally, we examine the scaling behavior predicted for 
the specific heat. In the clean, low-energy limit, 
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FIG. 4. Normalized vortex contribution to specific heat, 
SC e i(H)/[ lel Ta(H/H c2 ) 1/2 ] vs. Y = a(H/H c2 ) 1/2 T c /T for 
fixed temperatures T and scattering rate s F a s shown; unit of 
energy T c0 . Asymptotic large- Y limit \PlpK (dashed line). 

this result may be derived by substituting Eq. (4) 
into (5). One finds SC(H)/[j n Ta(H/H c2 )^ 2 } = F C [Y], 
where Y = a{H/ H c2 ) 1/2 T C /T and F c is a scaling func- 
tion which varies as F c ~ 31og2A V/(47rT c ) for Y < 1 
and as Fc — for Y 1. The full numerically de- 

termined, low-energy clean limit scaling function is plot- 
ted in Figure 4, along with full numerical evaluations 
of (5), which includes corrections beyond the low-energy 
approximation, as well as impurity effects. Scaling is ex- 
pected for a given data set provided H, T are such that 
Eh and T are both larger than the impurity scale 7. 
For example, in the clean case chosen (open symbols), 
scaling is obtained over the full range of Y, whereas for 
the dirty system (filled symbols) scaling has broken down 
completely. Junod |n| has recently reported data which 
crudely follow scaling predictions, but scatter in the data 
is too great to ascertain if the deviations follow the pat- 
tern predicted here. 

Conclusions. We have placed the elegant scaling argu- 
ments of Volovik |^|,|l8| and Simon and Lee ]3| regarding 
the specific heat of a <i-wave superconductor in magnetic 
field on a concrete foundation, introducing a simple for- 
malism capable of including both the effects of disorder 
and of energies comparable to the gap scale. We have 
shown that the density of states of a dirty d-wave sys- 
tem varies as H log H rather than the expected for 
the clean system. As a result, the predicted scaling of 
the vortex specific heat SC (H) with \/H /T breaks down 
in a well-defined fashion. As sample variability in the 
cuprates is notorious, it will be important to understand 
these deviations. We are in the process of studying sim- 
ilar effects in transport properties in a magnetic field. 

Note added: After submission of this paper we learned 
that Y. Barash et al. ]2(J had also obtained Eq. (3). 
We are grateful to G.E. Volovik for pointing out this 
reference. 
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